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Abstract 

We study the properties of Wilson loops in three dimensional non-compact U{\) gauge theories 
with global abelian symmetries. We use duality in the continuum and on the lattice, to argue that 
close to the critical point between the Higgs and Coulomb phases, all correlators of the Wilson loops 
are periodic functions of the Wilson loop charge, Q. The period depends on the global symmetry 
of the theory, which determines the magnetic flux carried by the dual particles. For single flavour 
scalar electrodynamics, the emergent period is Q = 1. In the general case of N complex scalars 
with a [/(l)^ -1 global symmetry, the period is Q = N. We also give some arguments why this 
phenomenon does not generalize to theories with a full non-abelian SU (TV) symmetry, where no 
periodicity in Q is expected. Implications for lattice simulations, as well as for physical systems, 
such as easy plane antiferromagnets and disordered superfluids, are noted. 
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I. INTRODUCTION 



Three dimensional abelian gauge theories have been a subject of intense study ever since 
Polyakov's demonstration of monopole induced confinement in compact electrodynamics[l|. 
The principal tool used to study confinement properties is the Wilson loop, which corre- 
sponds to an insertion of an external charge-anticharge pair into the theory. Area law for 
the Wilson loop indicates linear potential between distant charges and confinement. 

Non-compact three dimensional U(l) gauge theories with matter fields exhibit a different 
interesting physical phenomenon: existence of conformally invariant critical points. The 
prime example of such a theory is the single flavour scalar electrodynamics, the so-called 
Abelian Higgs model. In a certain region of parameter space, this theory exhibits a second 
order phase transition, which can be understood as being due to spontaneous breaking 
of the topological £7(1) $ global flux symmetry^, 0, 0, 0]. The order parameter for this 
symmetry is the monopole operator V(x), which creates a Nielsen-Oleson vortex with flux 
2tt. In the Higgs phase of the theory vortex excitations have a finite mass and the U(l)$ 
flux symmetry is unbroken. On the other hand, in the Coulomb phase of the theory the 
flux symmetry is spontaneously broken, with the photon being the corresponding goldstone. 
The phase transition can thus be visualized as being due to proliferation of vortices. Since 
the Z7(l)$ symmetry is the only global symmetry broken as one crosses the critical point, 
one might suspect that the phase transition in the Abelian Higgs model is in the (inverted) 
XY universality class and can be described in terms of a dual local theory of a dynamical 
vortex field V(x). This hypothesis is supported bvan exact duality between certain lattice 
versions of the Abelian Higgs and XY theoriesjl, Ijj]. 

The duality has been used extensively to study observables in the Abelian Higgs model, 
such as correlation functions of monopole, as well as magnetic field, operators. However, the 
behaviour of Wilson loops near the phase transition has largely escaped theoretical attention. 
This is not surprising: in the non-compact U (1) theory there is no linear confinement between 
external charges, so the prime motivation for studying the behavior of Wilson loops is gone. 
Nevertheless, as we shall show below, absence of confinement does not preclude interesting 
behaviour of the Wilson loop across the phase transition. 

In this paper, we discuss how to incorporate Wilson loops into the dual theory of the 
Abelian Higgs model using both symmetry arguments in the continuum and explicit duality 
transformation on the lattice. We find that a Wilson loop of charge Q in the direct picture 
gets mapped into an infinitely thin external flux tube carrying flux 2ttQ in the dual picture 1 . 
Since flux 2ir is invisible, one immediately concludes that the universal physics near the 
phase transition is periodic in charge Q of the Wilson line, with period Q = 1. This means 
that the behaviour of integer Wilson loops across the phase transition is non-universal: the 
length over which integer external charges are screened does not diverge as one approaches 



1 This fact has been previously noted in Rcf. [6| , using an argument slightly different from the one presented 
here; however, the consequences for critical properties of the Abelian Higgs model were not discussed. 
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the critical point. We note that the periodicity in the charge Q is emergent rather than 
fundamental: short distance physics of scalar QED is certainly not periodic in Q. 

In the second part of this paper, we generalize our discussion to what we shall call 
the "planar" theory: scalar QED with N flavours and a U{1) N ~ 1 global symmetry under 
independent phase rotations of scalar fields, as well as a symmetry under permutations of 
the flavours. This theory has N types of global vortices, which carry a fractional magnetic 
flux 2vr/AT0. Near its critical point, the theory is believed to be dual to a theory with 
N — 1 U(l) gauge fields and N flavours of scalars, which represent global vortices of the 
direct modelja, 0, 0]. We show that due to the fractional charge of the vortices under the 
flux symmetry, a Wilson line of charge Q in the direct theory gets mapped to an external 
magnetic flux tube with flux 2nQ/N in the dual theory. Thus, near the critical point, the 
physics of the N flavour model is periodic in the charge Q of the Wilson line with period 
Q = N. 

The emergent periodicity in the charge Q can be explicitly tested by lattice simulations. 
We suggest one observable for lattice simulations: the electric field produced by a straight 
temporal Wilson line at the phase transition. Like all universal observables, the electric 
field will be periodic in the charge Q of the Wilson line. We also present some explicit semi- 
quantitative predictions for the coefficient of the electric field based on the 1/M expansion 
of the dual theory. 

Finally, we would like to understand whether the phenomena described above generalize 
to U(l) gauge theory with N scalar flavours and a full SU(N) global symmetry. We shall 
argue that the answer to this question is no: the physics in the SU(N) symmetric model 
is not periodic in the charge Q of the Wilson line. In contrast to the situation in the 
theory with abelian global symmetry, the candidate dual degrees of freedom in the SU(N) 
symmetric model cannot be associated with local fields charged under the flux symmetry. 
As a consequence, no periodicity in the charge Q emerges. 

We would like to note that besides being of general theoretical interest, the behaviour 
of Wilson loops across phase transitions in non-compact U(l) gauge theories is important 
for a number of physical problems. For instance, it is believed that the phase transition 
from an antiferromagnetic Neel state to a valence bond solid (VBS) state on a square lattice 



is described by iV = 2 scalar non-compact QED|12|. Il3l |. The model with SU(2) flavour 
symmetry corresponds to a spin-rotation invariant system, while the model with an abelian 
U(l) global symmetry describes, the so-called, easy-plane antiferromagnet. The later model 
is also believed to describe the phase transition from a superfluid state to a disordered state 
in a theory of lattice bosons[14]. A missing spin impurity in this class of models is represented 



by a Wilson line in the gauge-theory description 15| . We have applied the dual description 



of the Wilson loops presented in this paper to study impurity induced VBS susceptibility in 



easy-plane antiferromagnets in Ref.[16|. The predictions of [16| may be explicitly tested by 
lattice studies of phase transitions in antiferromagnets and possibly by STM experiments 
on cuprate compounds. 

This paper is organized as follows. In section [III we describe how to incorporate Wilson 
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loops into a dual description of the Abelian Higgs model. In section, II II I we generalize this 
description to a model with N flavours and a ^/(l)^ -1 global symmetry. A discussion of 
the curious periodicity in the charge Q of the Wilson loop appears in section IIVI Section 
IVl contrasts the behaviour of Wilson loops in theories with abelian and non-abelian global 
symmetry. Concluding remarks are presented in section IVII 

II. WILSON LOOPS IN THE ABELIAN HIGGS MODEL 
A. Duality and Wilson loops 

It is well known that in three space-time dimensions, near its critical point, non-compact 
N = 1 scalar electrodynamics is dual to a theory of a complex (pseudo)scalar field with a 
global U(l) symmetry 

aiaa. The Lagrangians of these two theories £1X6 clS follows, 

Lqed = + ^ ~ * A M* + ™ 2 \z\ 2 + §M 4 (2.1) 

Lxy = |^| 2 + m 2 |V| 2 + ||V| 4 (2.2) 

Here is a non-compact gauge field, F M = e^xd^Ax is the magnetic field, and z and V 
are complex one component fields. The duality is understood as being true for the range 
of parameters where Lq ED has a second order phase transition (which at weak coupling is 
believed to occur for g/e 2 sufficiently large). One way to understand the duality is by noting 
that the phase transition in scalar QED is driven by spontaneous breaking of flux symmetry 
C/(l)$, which is precisely the global symmetry of Lxy- The order parameter for breaking of 
the flux symmetry is the monopole operator V(x) - that is the dynamical field of Lxy- As 
we know, to each continuous symmetry there corresponds a conserved current. In the case 
of flux symmetry of QED, this pseudo-vector current is just the magnetic field F^, which is 
trivially conserved in the absence of monopoles, d^F^ = 0. Let's introduce an external field 
Hfj, that would couple to this current, 

5L QED = iH^ (2.3) 

Suppose we are calculating some correlation function with insertion of a string of monopole 
operators {V qi (xi}} of charge qi at points Xi. The gauge field A M in the path integral is then 
subject to the condition, d^F^ = J2i 2vrgj5(a; — Xi). Then under the transformation, 

H^H^ + d^a (2.4) 
Sqed -> S QED + i J dxd^aF^ = S QED - i J dxad^F^ = S QED - 2iri'Y^q i a(x i ) (2.5) 
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Hence, by introducing the field we can enlarge the global £7(1)$ symmetry to a fictitious 
local symmetry, provided that the monopole operators transform as, 



V q (x) -> e 2mqa{x) V q {x) (2.6) 

The dual Lagrangian Lxy has to posses this local symmetry. Hence, to introduce the field 
if M into the dual Lagrangian we simply have to covariantize the derivative of the dynamical 
monopole field V, 

d^V -> DpV = (0„ - 2mH li )V (2.7) 

in eq. (12.21) . Other "gauge invariant" operators can also be added to Lxy, e.g. H*; however, 
their contribution will, generally, either cancel out in correlation functions or be less singular 
near the critical point (see section Hi CI for a more detailed discussion). 

Thus, the dual Lagrangian in the presence of a background source field is given by, 

Lxy = |(<9 M - 2mH fl )V\ 2 + m 2 \V\ 2 + ||Vf (2.8) 



The covariantization procedure (12.71) was explicitly written down in Ref. [l7|. Similar 
arguments for the case of a constant imaginary H^, which physically represents an external 
magnetic field in the QED language and translates into a chemical potential for the flux 



symmetry in the XY language, have been given in Ref. [181. In the next section, we shall 
also give an argument based on an exact duality transformation on the lattice, which will 
support (12.81) . 

Having learned how to incorporate the source field into the dual Lagrangian, it is now 
trivial to dualize Wilson loops. Indeed, insertion of a Wilson loop W(C) into a correlation 
function is equivalent to adding into the Lagrangian the source term 

5L = iQ [ dx^ = iQ [ dS^ = i [ dxH^ (2.9) 
Jc Js J 

where the surface S satisfies, dS = C and 

H ll {x) = Q( dS,5{x~y) (2.10) 
Jyes 

So is a field that lives on the surface of the Wilson loop and is directed perpendicular to 
this surface. 

Another benefit of introducing the source field is that by differentiating with respect 
to it we can compute correlation functions of the magnetic field F^. For instance, 

HW)b = 6l sH^ x f = - 2 ™((V^D„V - (D^V) {x)) H (2.11) 
Hence the topological flux current F„ of QED gets mapped into the Noether's current asso- 
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ciated with the global U(l) symmetry of the dual model. 

We have seen above that the Wilson loop W(C) in the dual theory is specified by a surface 
S with dS = C rather than by the contour C alone. Let's investigate the dependence of the 
dual theory on the choice of this surface. If we pick a different surface S', with dS' = C then 
the field undergoes a gauge transformation — > H' = + d^a with a(x) = —Q l ie v 
where V is the volume bounded by the two surfaces S and S' . Hence, 

(V(x)...) HI = e 2 ^(V(x)...) H (2.12) 

where ellipses denote some other operators. Thus, the operator V(x) is invariant under 
changing the surface of the Wilson loop if and only if Q is an integer. This is nothing but 
Dirac's condition expressed in the language of the dual theory. However, a theory with 
arbitrary non-integer Q is still sensible provided that we don't consider monopole operator 
insertions, or more formally, confine our attention to correlation functions of operators which 
are invariant under the fictitious U(l)$ local symmetry, e.g. the magnetic field operator 
— iF^ = —l-niV^D^y. In fact, if we are dealing with such gauge invariant operators we 
don't necessarily have to use the precise form of H given by f 1 2 . 1 j) ; defining 7 M to be a field 
living on the perimeter of the Wilson loop and directed along it, 

1^x) = q! dy^x-y) (2.13) 
Jyec 

we see that, 

e^uxd u H x = 7m (2.14) 

Then, by performing a suitable gauge transformation on and V we can choose to be 
any field with curl given by 7 M . Thus, we see that the duality maps a Wilson loop of charge 
Q in the QED language to an external magnetic flux tube of flux 2wQ in the XY language. 
This correspondence has been noted in Ref. @, but the consequences of this correspondence 
for the critical properties of Wilson loops were not discussed. 

Thus, we have to solve an Aharonov-Bohm like problem for the dual field V. The question 
is simplest to analyze with the original gauge choice (I2.10p . This "string" gauge is equivalent 
to Hp = and the boundary condition, 

V(x + ) = e 2niQ V(x~) iorxeS (2.15) 

where x denote points on opposite sides of the surface S (x = x±en, for e — > + , where n is 
a local normal to S). So the Wilson loop imposes a twisted boundary condition (I2.15P in the 
dual theory. We observe that the physics is, therefore, a periodic function of Q. For integer 
Q the boundary condition ( 12. 15ft is trivial - there is no twist. So our argument indicates 
that integral Wilson lines do not affect the physics on distances of order of the correlation 
length of the theory: screening of integral charges takes place on length scale which does not 
diverge as one approaches the phase transition. This is certainly an unexpected result: we 
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will discuss it further in section IIVI However, first we would like to obtain further support 
for this result by performing an explicit duality on the lattice, which will provide additional 
physical insight into the origin of the periodicity in Q. 

Another interesting consequence of periodicity in Q is the emergence of charge conju- 
gation C and charge parity CP invariance at points Q = ±1/2. Indeed, the Wilson line 
generally breaks both C and CP symmetries, which map the charge of the Wilson line 
Q — ► —Q. However, due to periodicity in Q, the points Q — ±1/2 are identified, so C and 
CP symmetries are effectively restored for half-integer valued Q. 

B. Duality on the Lattice 

In Section III Al we have given arguments on how to perform the duality on QED 3 with 
Wilson loop insertions. Our arguments were very general, being based on the presence of 
flux-symmetry alone. In the present section, we would like to support the arguments of 
III Al by performing an exact duality between lattice versions of QED and XY model. In 
the process, we will obtain some insight into the emergent periodicity in the charge of the 
Wilson line. 

The lattice duality between non-compact QED and XY model is very well known [3, 0] . 
We start from the QED lattice action, 

S Q™ = h + Ya ^ {d6 (2 ' 16) 

~m it* 

Here A, M is a gauge field living on links of the lattice and e td is a matter field, whose amplitude 
is frozen. The auxiliary variables integers living on the links of the lattice, whose 

purpose is to ensure the 2n periodicity of the variable 9. For our purposes it will also be 
useful to add sources corresponding to monopole and Wilson loop insertions into the action, 

S QEB = ^2 E( QA - 27rr )l + kg ^ {d0 ~ A ~ 2?m) ^ + 1 ? H '^ UA ~ 27Cr)] » (2 - 1?) 
m w m 

Here rj M is an integer valued source field representing Dirac strings running from locations 
of monopoles to infinity, satisfying, 

i 

where and ji are correspondingly charges and locations of monopole insertions. The field 
in (I2.17P is a source coupling to the physical magnetic field (OA — 27ir)j fl (i.e. the 
magnetic field with the Dirac string subtracted). To represent a Wilson loop of charge Q, 
we can choose Hj„ to be equal to Q on the surface perpendicular to the loop and zero 
everywhere else. For integer Q, such a Wilson loop will be independent of the choice of 
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the surface, while for non-integer Q it will depend on the choice of the surface if monopole 
operator insertions are present. 

Now, we perform the duality. First, we decouple the kinetic term for the gauge field by 
introducing an auxiliary field Pj M . We also Poisson resum the field n JAt by introducing an 
integer valued variable Jj- M . 

w-yE4+^E(^- A - + 2m E n ^ + 1 E( p + H )U nA - 2 ^)m 

3H J> J> in 

(2.19) 

After Poisson resummation, n« a becomes a free real variable, so we can shift it, 27m' = 
(2ixn + A - dff)^, 

~m w it* ~m 

(2.20) 

Now, performing the integral over we obtain a constraint 

V-J = (2.21) 

Physically, Jj M represents the world-lines of e ld particles. If we were to integrate over all 
other fields in the problem, we would see that these worldlines interact with a long-range 1/r 
Coulomb interaction (there are also local interactions between these worldlines controlled 
by the coupling strength g). We solve the constraint (12.211) in terms of an integer valued 
field b- jfl , 

J m = {Ub) m (2.22) 
Now, we rearrange our action slightly and integrate over the n' field, 

Sqed -yE P i + 8 E ^ P + H ) ~ J h» A ^ ~ 2m E( P + H ^ + f E 4 ( 2 - 23 ) 

ip- in it* 

Performing the integral over the A field, we obtain a constraint, 

□(P + H) jfM - J jfl = (2.24) 
Recalling (12.22ft we can solve ( 12. 24ft by introducing a real field (p, 



P + H -b = ^ (2.25) 
2n 



arriving at the action, 



2 

Sqed - ^ J>y> ~ 2nH + 2nb ^l + f E( D C ~ 1 E^ + 27r %^> ( 2 - 26 ) 
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Recalling that both b and r are integer valued, summing the last term by parts and using 

(GOD, 

Eq. (12.271) is the final form of the dual lattice action. Temporarily setting g = 0, dropping 
all the source field, we obtain the usual Villain form of the XY model, with e lip being the XY 
field. As already noted, the role of finite g is to introduce short-range interactions between 
e particles, that is between vortices of the ip field. This is also evident from the dual action 
(12.271) as we can identify Ob with the density of ip vortices. Since such vortices for e 2 ^ 
already interact with a long-range Coulomb potential, we expect the phase transition at 
finite g to be in the same (inverted) XY universality class as at g = 0. 

Now, restoring the source fields into (I2.27p . we immediately identify e lip , with the 
monopole field of QED. Moreover, we also see that the source field H enters the dual action 
by gauging the lattice derivative of the ip field. We have predicted this fact from symmetry 
arguments alone in Section III Al Hence, we see that the action ( 12.81) is a suitable continuum 
generalization of our dual action (I2.27p . where we identify V ~ e l{f . 

Now, we come to the question that interests us most: what is the influence of integral 
Wilson loops on our theory. For a general non-integral charge Q the Wilson loop enters the 
dual theory as a highly non-local coupling, as the source field H lives on the whole surface 
of the Wilson loop, rather than on its perimeter. Nevertheless, for an integral charge Q, we 
can perform a transformation, 

b' = b-H (2.28) 
as in this case H is integer valued. Hence, dropping the monopole insertions, 



6 - Y^dtp + 2nb% + | + !)% (2.30) 



where we used (□if)j A1 = 7j M with the Q valued vector field 7 pointing along the perimeter 
of the Wilson loop (and being zero everywhere else). Thus, for integral Q the coupling of 
the theory to the Wilson loop becomes local. So we expect the physics of integer Wilson 
loops to be drastically different from that of non-integer ones. 

Moreover, we see that for the special value g = 0, the coupling to the Wilson loop in 
( 12.30!) disappears all together. From the point of view of the direct theory ( 12.161) this fact 
is not surprising: setting g = forces A = d9 — 2irn, so that W(C) = 1. Thus, the theory 
possesses a limit in which the insertion of integer Wilson loops is trivial, but as one varies 
e 2 , a phase transition in the (inverted) XY universality class still occurs. As already noted, 
we expect the phase transition at g 7^ to be in the same universality class as that of g = 0, 
that is g is in some sense an "irrelevant" coupling. Hence, one can argue that the physics of 
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integer Wilson loops is non-universal - their only effects appear on distance scale of order of 
the microscopic cutoff of the theory. 

To illuminate this conclusion further, let us rewrite the dual theory (I2.27P in terms of 
magnetic vortex lines, (i.e. worldlines of the field e %lf ). We proceed by Poisson resumming 
the variable bj^ by introducing an integer valued field Ij^, 

s q* d = ^5 E(^ - 2nH + 2n ^h + 1 E( Q C + 2m E 'jA. (2- 31 ) 

where we've dropped the monopole source s for simplicity. Now, shifting 2-Kb 1 = 2-n;b — 2iiH + 
dtp, 

Sqed - y X) 6 ^ + 1 E( Q6 ' + ^> + < E W 2 * 6 ' + M - <H> ( 2 - 32 ) 

Integrating over the field we obtain the constraint, 

V-/ = (2.33) 

Physically, are just the magnetic flux tubes. We can also integrate over the field b', 
obtaining, 

Sqed = \ Yp^Y^iA^Y^ + 2 ™ E H 3^ D n')h, + ^ E ^^"7,> (2.34) 
ii'M w'a« jj'm 

where the propagator Djji is given by, 

D =Ly 1 ^_ e «Cf-iO (2.35) 

with a being the lattice spacing and V^-the number of sites in the lattice. The last term in 
(I2.34p does not couple to the dynamical field I and, thus, is trivial. The first term in (I2.34p 
represents the short range interaction between the flux-tubes. The second term, 

z$ = 2m X H- m {e 2 D n ,)l JV (2.36) 

is the one that interests us in conjunction with the properties of Wilson loops. Indeed, we 
identify, 

B ] , = 2^(e 2 %^^ (2-37) 

r 

with the magnetic field produced by each flux tube. We see that the magnetic field due to 
each flux-line is short-range, as expected. Hence, if we take H to represent the Wilson loop, 
(I2.36P simply adds up the contribution of all flux tubes to the flux through the loop. In the 
special limit g = 0, the propagator Djji is ultra-local, e 2 Djji = Sjji and the flux-lines become 
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infinitely thin. Then, (12.361) simply counts the number n of flux-lines passing through the 
loop (in other words, n is the linking number of the flux-tubes with the Wilson loop), 

z$ -> 2mQn (2.38) 

We immediately see that for Q - integer, the term (12.381) gives a trivial contribution. 

Now, turning g back on, our flux-lines obtain a finite thickness r/. Then the expression 
for the flux (12.381) is modified by contributions from flux-tubes passing within a distance 
~ rf from the boundary of the Wilson loop (that is flux-lines, which are not entirely inside 
or outside the loop). This contribution can be understood as a local coupling to the Wilson 
line, which is expected to be less relevant in the RG sense than the non-local term (12.381) . 
In fact, in the next section, we will argue that all such local, line-like coupling terms are 
irrelevant. 



C. Perturbations local at the Wilson line 

Another way to argue the non-universality of response to integral external charges is to try 
to construct a relevant perturbation of the continuum dual theory (12.81) . Such perturbations 
have to be invariant under the fictitious £7(1)$ local symmetry discussed in Section III Al 
Moreover, as discussed in section Hi B\ for integral Q, the perturbation must be local to the 
Wilson line. There are plenty of operators with correct symmetry properties, since in the 
action (12.81) we included only the most relevant (in terms of power counting) terms. We 
count the source field as having dimension 1. This is the canonical dimension of this 
field. Taking into account C and P symmetry, the operators of lowest dimension which we 
can write down are, 

[{t, vX d u H x ) 2 ] = [ 7 g = 4 (2.39) 
[{e^ x d v H x ) W] = [llV^V] = D + 2 (2.40) 
[{-ie»uxD v V^D x V){e ml} d a Hp)\ = [{-ie^D^ D x V) lfl ] = D + 2 (2.41) 

where we've included the canonical dimensions of the operators in question. The first op- 
erator (12.391) does not couple to the dynamical fields of the dual theory and, therefore, its 
contribution is trivial. The other operators are irrelevant by power-counting. We particu- 
larly want to draw attention to the operator (12.41 j) . since it couples the Wilson line 7m to 
the vortex density operator, — ie^DyV^ D X V . This is precisely the kind of coupling that 
we have for integral charges in the dual lattice action (12.301) . namely Ob >r y. 

We note that the above argument might be too naive, as the field strength 7ai = 
£[iu\d u H\ ~ S 2 (x) corresponding to the Wilson line is very singular. For instance, it is 
not clear how to interpret the j 2 term in (I2.40p . If one takes 7^ ~ (5 2 (x)) 2 ~ A 2 5 2 (x), then 
the extra factor of A 2 upsets our power counting scheme. Such singularities will occur at all 
orders in field-strength. So perhaps it is more appropriate for integer-valued external charge 
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to simply write all operators in the dual theory, which live locally on the Wilson line (of 
course, we then don't know the dependence of the coefficients of these operators on charge 
Q). The operator which would seem to be most relevant is, 



5L 



u 



drV^V(x = 0,T) 



(2.42) 



Here u is some real coupling constant, and at the critical point its scaling dimension is 

dim[w] = 1 - dim[V^V] = 1 - (3 - l/v) = 1/u - 2. Because 



determined following Ref. [19 
v ~ 2/3 for the XY model, u is an irrelevant perturbation. We, therefore, come to the 
curious conclusion that it is impossible to construct a (weak) perturbation localized on a 
one-dimensional line, which would be relevant at the fixed point of the XY model. 



III. WILSON LOOPS IN Y-FLAVOUR MODEL WITH U(l) N ~ l SYMMETRY 

In this section, we consider a theory with N flavours of scalar fields z a (N does not 
necessarily have to be large), 

L = ^1 + - iA,)z a \ 2 + U{z a ) (3.1) 

We refer to this theory as the "planar" model. Here, U is some potential with the global 
U(1) N symmetry under independent phase rotations of the z a fields. The singlet component 
of this symmetry is actually gauged by the field A^, 

U(l) : z a - e w ^z a , A,^A, + d,6 (3.2) 

while the non-singlet components are true global symmetries of the theory, 

Uilf- 1 : z a - e ieata "z a (3.3) 

where t a , a = 1..N — 1 are the generators of the [/(l)^ 1 symmetry satisfying, Y2 a ^a = 0- 
We require U to have a symmetry under the permutation of labels of z a fields. We choose U 
in such a fashion that in the "condensed" phase of the theory, it favours non-zero expectation 
values of all components of the z a field, so that the vacuum manifold of the theory is a torus, 
(S 1 ) 1 ^ (here we temporarily forget that the singlet symmetry is gauged). For N = 2 the 
theory under consideration is believed to describe the phase transition in the easy-plane 
ant if erromagnet . 

We would like to dualize the theory (13.11) . Similar theories were dualized in Refs. 
and here we will present a related discussion aimed at incorporating Wilson loops into the 
dual theory. An exact duality on the lattice appears in the appendix, but we can write down 
the form of the dual action from very general considerations. Let us first identify the dual 
degrees of freedom. We go to the condensed phase of the theory (13. ip . where all (z a ) ^ 0. 
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Then, we can have vortices in any component of the z a field. Formally, the homotopy group, 
ni((S l ) N ) = 1i N . So, we have iV types of vortices, which become the degrees of freedom of 
the dual theory V a , a = 1..N. 

These vortices are global, rather than local. Indeed, let's consider a vortex in the first 
component zi, 

z 1 (x) ~ ve tX ^\ z a ~v,a^l, \x\ — > oo (3.4) 

where X(x) winds from to 2tt as one goes around a contour out at infinity surrounding 
the vortex. Then, this vortex corresponds to a space-time dependent transformation of the 
vacuum 1EE2|) . (l3~3l) . with, 9{x) = ±\(x) and 6 a {x)t a = (1 - 1/N, -1/N, ... - l/N)X(x). 
Thus, our vortex possesses a winding both in the local and in the global symmetry group. 
The winding in the local U(l) group will be canceled by the gauge field, 

M*) = 9 A*) = ^X(x) (3.5) 

hence our global vortices carry a magnetic flux $ = 2tt/N. Therefore, under the flux 
symmetry (12.41) . the fields V a should transform as, 

Vjx) -> e 2ma ^ /N V a (x) (3.6) 

This fact will be crucial for the analysis to follow. 

The winding in the global group will lead to a long-range Coulombic interaction between 
our vortices. We will need dynamical gauge fields in the dual theory to give rise to this 
interaction. However, if we have a unit winding in each component of the z field, our vortex 
becomes completely local, and carries total flux 2ir. We can think of such a local vortex as 
a composite of N global vortices of different types. The creation operator for this flux-tube, 
therefore, will be, 

V(x)=l[V a (x) (3.7) 

a 

Since the local vortex carries flux 2tt, we can also associate the operator ( 13. 7\\ with the 
monopole operator of the direct theory. Indeed, given (13.61) . under the flux symmetry ( 12.41) . 

V(x) -> e 27Tia{x) V(x) (3.8) 

which is the correct transformation law for the monopole operator (12. 6p . 

We expect local vortices to interact by short range forces. Therefore, the operator ( 13. 71) 
should not be charged under the emergent gauge fields of the dual theory. 

We are now ready to write down the dual theory, 

1 2-7T? 
i 

Here = B a J, a a , a = 1..N—1, are emergent dual gauge fields, which couple to the non-singlet 
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currents. F a = e^xdyB" are the corresponding field strengths. The dual potential U(V a ) 
is chosen to have the same properties as the direct potential U: it has a U(1) N symmetry 
under independent phase rotations of the fields V a and a symmetry under permutation of 
labels of V a fields. Moreover, it favours (V a ) ^ for all a in the condensed phase of the 
dual theory. Thus, the theory (13.91) has a local {/(l)-^ -1 symmetry, 

Uilf- 1 : V a (x) ^ e l ^ x ^V a (x), B% ^ B% + d^ a (3.10) 

as well as the global U(l) flux symmetry of the direct theory (13.61) (which we have promoted 
to a local symmetry by introducing a non-dynamical source field H^). As required, the 
monopole operator (13.71) is invariant under the local [/(l)^ -1 symmetry of the dual theory 

The theory (13.91) also has a global U(l) N ~ 1 symmetry associated with conservation of 
fluxes of the N — 1 emergent gauge fields. This topological symmetry can be identified with 
the Noether's symmetry (13.31) of the direct theory. 

Now, we would like to apply the duality discussed above to study the properties of Wilson 
loops. Recall, that to represent Wilson loops we must use a source field given by (12.101) . 
As discussed for the case of iV = 1 theory, the effect of such a source field on the dual action 
(13. 9p is to introduce a twisted boundary condition for the vortex fields, 

V a {x + ) = e 2wiQ/N V a {x-), forxeS (3.11) 

where Q is the charge of our Wilson line. The physical origin of the factor 1/N is the 
fractional charge 2n/N of the vortex fields V a under the flux symmetry. Thus, we come 
to the amazing conclusion that the universal physics in the planar model is periodic in the 
charge Q of the Wilson line, with period Q = N. This is a generalization of the Q = 1 
periodicity of single flavour QED discussed before. 

Similarly to the N = 1 case, a consequence of the Q mod N periodicity is the emergence 
of the C and CP symmetries for half-integer values of Q/N. As before, this periodicity is 
due to the identification of C, CP conjugate points Q/N = ±1/2. 

IV. DISCUSSION 

We have argued above that non-compact QED with N identical flavours and a {/(l)-^ -1 
global symmetry acquires a Q mod N periodicity in the charge Q of the Wilson loop near its 
critical point. In the particular case of iV = 1, Abelian Higgs model, the period is Q = 1. It 
is important to note that this periodicity is emergent, rather than fundamental. On shortest 
distance scales, QED in three dimensions is perturbative; the electric field produced by an 
external charge Q is simply Coulombic, E = which is obviously not periodic in Q. 

One also should not confuse the emergent periodicity with trivial screening of exter- 
nal charges by dynamical fields. Such screening generically takes place in low-dimensional 
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abelian gauge theories and is due to the confining nature of the Coulomb potential (linear in 
one spatial dimension and logarithmic in two), which leads to binding of a dynamical particle 
by the external charge 2 . However, this trivial screening typically i) occurs on distance scales 
r ^> £, where £ = m -1 is the correlation length of the theory, ii) leads to a resulting period 
of Q = 1. One classic example of this phenomenon is the Q mod 1 periodicity of string 



tension in massive Schwinger model[20|. The phenomenon considered in the present paper 
is clearly different since i) the screening occurs on distance scale r £ (for strong coupling, 
e 2 ~ A, the screening is actually expected to take place on short-distance cutoff scale), ii) 
the resulting period is Q = N rather than Q = 1. For N > 1, in the Coulomb phase of 
the {/(l)-^ -1 symmetric theory, we also expect the usual screening on distances r ^> £ with 
period Q = 1. 

Physically, the periodicity discussed in this paper is due to the fact that the degrees of 
freedom responsible for the phase-transition are local fields carrying a fixed magnetic flux, 
$o = 27r/N. The Wilson loop of charge Q is then expressed through the linking number n 
of the worldlines of dual particles with the contour of the loop, 

W(C) = e iQ *° n (4.1) 

and is trivial for Q = mod N. For the special case of N = 1 this phenomenon is a 
manifestation of the fact that the phase transition is driven by magnetic rather than electric 
degrees of freedom. 

The results presented in this paper can be explicitly checked by lattice simulations. The 
simplest lattice counterparts of the continuum theories under consideration were actually 
discussed in this paper in sections III Bl and in the appendix. These lattice theories were 
previously simulated in a number of studies^, @, 11, 17]. We have predicted that near 



the phase transition all physical observables become periodic in Q. However, we have not 
discussed specific observables. One observable that we suggest for lattice simulations is 
the electric field produced by a straight temporal Wilson line (we define the electric field 
Ei = F i3 = —EijFj). At the critical point, the electric field must have the form, 

(-z£7,) = C(Q)l (4.2) 

This form is dictated by the fact that the magnetic field — iF^ is a conserved current, which 
receives no renormalization and hence has conformal dimension D — 1 = 2. The coefficient 
C(Q) is a universal function depending only on the charge Q of the Wilson loop 3 . We predict 
that C(Q) must be periodic with period Q = N. In particular, C(Q) vanishes for Q = 
mod N. Moreover, due to the emergence of C and CP symmetries at half-integer Q/N, 



2 In two dimensional theories, an infinitely large Wilson loop of charge Q is equivalent to a non-zero value 
of the topological 9 angle, 9 = 2ttQ. So the large-distance periodicity in Q is synonymous to periodicity 



m 

3 



C(Q) is a real valued function: the electric field is imaginary as we are working in Euclidean space. 
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FIG. 1: Coefficient C(Q) of the electric field at the critical point of the Abelian Higgs model, see 
eq. (I4.2|) . computed from a 1/M expansion of the dual theory. 



C(Q) actually vanishes for both integer and half-integer Q/N. This does not mean that 
the electric field vanishes at these special points. Rather it will be controlled by irrelevant 
couplings, e.g. one in eq. (12.411) . and will fall off as some higher power of 1/r with a 
non- universal coefficient. 

The present paper has concentrated on incorporating Wilson loops into the dual La- 
grangian and discussing their general properties. Nevertheless, the dual Lagrangian can also 
be used for explicit calculations of properties of Wilson loops. Of course, the dual theory is 
still strongly coupled in the infra-red: we have mapped one difficult problem onto another. 
However, in the case of the Abelian Higgs model (N = 1), the dual theory is just the three 
dimensional U(l) symmetric scalar field theory rather than a gauge theory. A wealth of nu- 
merical and analytical information is known about the phase transition in this theory (XY 
universality class). It is known that e-expansion and large-M expansion (whereby the dual 
field V is promoted to have M components), produce accurate results in this theory. This is 
in contrast to the e expansion in the direct theory, which predicts the existence of a critical 
point only for N > 182 (see Ref. 21J) (here the field z is promoted to an N component 
SU(N) multiplet). Likewise the large- N expansion of the direct theory when extrapolated 
to A" = 1 is known to produce results for anomalous dimensions, which are numerically 
notoriously inaccurate. Moreover, we shall argue below that the SU(N) symmetric theory 
actually does not possess any periodicity in the charge Q of the Wilson line and hence does 
not capture the qualitative features of the N = 1 theory. Thus, at least for the N = 1 case, 
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there are clear advantages of performing calculations in the dual, rather than direct theory. 
In Ref. Ila, we have used the large-M expansion of the dual theory to explicitly compute the 
universal function C(Q) of eq. (14. 2p at M — oo, see Fig. [TJ The coefficient A of the linear 
term of the expansion of C(Q) ~ AQ for Q — > is actually related to the conductivity in 
the XY model and is known numerically from Monte-Carlo simulations A pa 0.29, as well as 
to Q(e 2 ) in e expansion A pa 0.32, and to next to leading order in 1/M expansion A pa 0.25 
(see |l6j and references therein). 

For {/(l)^ -1 symmetric theories with N > 1 the dual theory is a gauge theory with 
iV — 1 gauge fields (the N = 2 theory is actually self-dual). Thus, explicit calculations in the 
dual theory are unlikely to be numerically accurate. However, they may illuminate general 
features, which are not immediately obvious in the direct theory. 



V. WILSON LOOPS IN THE SU(N) SYMMETRIC THEORY 

It is interesting to ask whether the periodicity in the charge Q of the Wilson loop gener- 
alizes to theories with N flavours and a full SU(N) global symmetry. We shall argue below 
that the answer to this question is no. 

A powerful tool in the analysis of SU(N) symmetric model is the 1/N expansion, which 
allows one to study the physics directly without performing any duality transformations. 
The 1/N expansion is typically performed in the limit e 2 — > oo, so that the bare kinetic 
term for the gauge field is absent. One also usually replaces the short range repulsive 
interaction between scalar fields by a hard-constraint, Y2 a z a z a — ~> which can be enforced 
by a local Lagrange multiplier A. One then obtains the Lagrangian of the CP N ~ X model, 

L=\{d,-iA^)z\ 2 + i\{\z\ 2 --) (5.1) 

9 

The behaviour of Wilson loops with Q ~ 0(1) in N is easily captured by the 1/N expansion. 
For instance, at the critical point of the theory one recovers the form of the electric field 
(14. 2p . where to leading order in 1/N, C(Q) = The SU(N) symmetric theory, thus, 
clearly does not display a Q mod 1 periodicity of the single flavour theory. One might not 
be too surprised by this fact, since we already saw that the iV-flavour theory with &U(1) N ~ 1 
symmetry has only a Q mod N periodicity. Does the SU (N) symmetric theory share this 
periodicity of its planar counterpart? 

To answer this question we must take Q ~ O(N) in the large N expansion. In this limit, 
the Wilson loop will modify the saddle-point of the expansion. For iV = oo, at the critical 
point, we should no longer expand around A^ = 0, A = 0, but rather around space-time 
dependent A^ and A. By dimensional analysis, a straight temporal Wilson line placed at 
the origin will generate saddle point fields, 

M,(r) = fA i\(r) = b M (5.2) 
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The coefficients a(Q) and b(Q) should be chosen in such a way that saddle-point equations 
are satisfied, 



(J^x)) = (z^^x)) = -Q8^8 2 {x) (5.3) 
(z*z(x)) = - (5.4) 

The expectation values in (15.31) . (15. 4p are to be computed by finding the propagator of the z 
fields in the background of the saddle point fields (15.21) . We have not been able to find this 
propagator explicitly. Nevertheless, the saddle-point equation ( 15. 3ft is not periodic in Q/N, 
thus, we conclude that the SU(N) symmetric theory is not periodic in Q. 

One may ask, what makes the SU(N) symmetric theory so different from its deformation 
with a £7(1) 1 symmetry. Our results regarding the planar theory relied on rewriting 
the problem in terms of dual degrees of freedom. Such a duality transformation, either 
in the continuum, or on the lattice is, so far, not known in the SU(N) symmetric theory. 
Nevertheless, we would like to provide some speculations on the possible form of the dual 
theory and implications for the properties of Wilson loops. 

To construct the dual theory, we first need to identify the dual degrees of freedom. For 
the planar theories considered above these are vortices: pointlike topological defects of the 
two dimensional reduction of the theory. Then the three dimensional theory is formulated 
in terms of vortex loops, i.e. worldlines of the two-dimensional defects. We would like to 
follow the same procedure for the SU (N) symmetric theory: as a first step we need to find 
two dimensional topological defects. Fortunately, classical defects of the two dimensional 



CP 1 model are the very- well known instantons [22]. These instantons carry magnetic 
flux $ = 2nq, where q is the integer topological charge of the instanton. To make further 
connection with the planar iV-flavour model we recall that all the solutions with topological 
charge q are known exactly, and can be parameterized in terms of Nq complex coordinates, 
a a i, with i = l..q. For q > 0, 

<? 

w a (s) = c a Y\_(s-a ai ) (5.5) 
i=i 

**(*) = , ( 5 - 6 ) 

{W'W{S))2 

where s = x\ + ix<i- The coefficients c a specify the overall orientation of the instanton in 
flavour space at infinity and are not very important. The corresponding expressions for 
q < can be obtained by taking s — > s. We see that the variables a a i are locations of 
vortices of fields z a . Thus, just as in the case of the planar iV-flavour model, the instanton 
(vortex) with flux 2-irq can be decomposed into qN fractional instantons (vortices) to which 
we can assign flux 2n/N. 

However, there is one major distinction between the planar and SU(N) symmetric mod- 
els. For the planar model the magnetic flux is concentrated in the core of the fractional 
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vortices, which are assumed to have some microscopic size. On the other hand, for the 
SU (N) symmetric model, there is no notion of the core-size since the theory is (classically) 
conformally invariant. The flux density of a configuration of fractional instantons is not 
localized near their positions dj a , but rather is smeared out in a distribution that depends 
in some highly complicated manner on the ratios of distances between these constituents. 
An instructive example is the instanton with charge q = 1, where the flux produced by the 
fractional instantons always clumps together into a rotationally invariant distribution. 

Thus, despite their similarity to fractional vortices of the planar model, fractional instan- 
tons do not carry a local flux and upon transition to three dimensions cannot be promoted 
to local fields charged under the £7(1)$ symmetry. The expression for the Wilson loop (14. ip 
in terms of the linking number of vortices with the loop contour is, therefore, inapplicable for 
the SU (N) symmetric case, and no periodicity in the charge Q of the Wilson loop appears. 



VI. CONCLUSION 

The purpose of this paper was to incorporate Wilson loops into the dual description 
of critical non-compact abelian gauge theories in three dimensions. This goal has been 
achieved for non-compact QED with iV flavours of identical scalar fields and a [/(l)^ 1 
global symmetry. The Abelian Higgs model corresponds to the N = 1 case of our general 
construction. A remarkable property, which follows from the dual description, is that the 
universal physics close to the phase transition is periodic in charge Q of the Wilson loop 
with period Q = N. In section [TV] we have provided a detailed discussion of this unexpected 
result; here we repeat a few conclusions. We have argued that the periodicity is emergent 
at the phase transition rather than fundamental to the theory. We also claim that this 
behaviour is distinct from trivial screening of electric charge in the Coulomb phase of low- 
dimensional abelian gauge theories. Moreover, we've argued that this periodicity does not 
generalize to the theory with a full SU(N) invariance. Thus, any attempt to understand the 
behaviour of Wilson loops at the phase transition of the Abelian Higgs model through the 
1/N expansion of its iV-flavour, SU(N) symmetric counterpart will fail to reproduce this 
qualitative feature. 

The predictions of the present paper can be explicitly tested by lattice simulations. In 
section [TV] we have suggested one observable: the electric field produced by a straight 
temporal Wilson line at the critical point, to test the periodicity discussed above on the 
lattice. 

Finally, the application of results of the present paper to the description of an impurity 
in a two-dimensional antiferromagnet in the neighbourhood of the quantum phase transition 



to a valence bond solid (VBS) state is a subject of a separate work [16|. It has been argued 



in 



1 61 ] that a missing spin impurity induces a vortex in the VBS order parameter. We hope 



that such vortices will be observed by future STM experiments on cuprate compounds. 



19 



Acknowledgments 



The present work evolved out of collaboration with Subir Sachdev aimed at understanding 
impurities in quantum antiferromagnets 16j. I am grateful to Subir for many useful discus- 
sions and comments. I would also like to thank Ariel Zhitnitsky for sharing his expertize on 
fractional instantons. This work was supported by NSF Grant No. DMR-0537077. 



APPENDIX A: DUALITY IN THE PLANAR THEORY ON THE LATTICE 

In this appendix, we discuss a lattice counterpart of the conjectured continuum duality 
discussed in section IIHI We start from the action, 

s q™ = ^Y.( nA - 2 ^l+^ g Y.Y.( dda - A - 2 ™ a ^ ( A1 ) 

This is a generalization of U(l) lattice theory (12.171) . in which we have introduced N flavours 
of dynamical matter fields e l6a , with corresponding integer valued variables n a that ensure 
the periodicity of 8 a variables. As before, we decouple the kinetic term for the gauge field 
by introducing an auxiliary field Pj^ and we Poisson resum the fields by introducing 
integer valued variables J? , 



(A2) 

Shifting, 2vm ,a = 2?™° +A-d6 a , 

jfi a j/j, a j/j, jn 

(A3) 

Performing the integral over 9" we obtain a set of constraints 

V • J a = (A4) 

Physically, Jj* M represent the world-lines of e 10 " particles. The singlet combination J2 a J a nas 
long-range 1/r Coulomb interaction, while the non-singlet combinations have short range 
interactions. We solve the constraints ( ]A4j) in terms of integer valued fields 6? , 

J« = (□&■%> (A5) 
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Integrating over the n' fields 



v^2S p J +l S( D ( p +^-E n,A> - 2m E( p + + ~ E E J ^ 2 

(A6) 

Performing the integral over the A field, we obtain a constraint, 

□(P + H)j, - = (A7) 

a 

Recalling (IA5I) we can solve (I AT[) by introducing a real field (p, 

a 

arriving at the action, 

s - £-2 - n> + § E B D6 % - 4 E ^ ( A9 ) 

We now Poisson resum the field b a by introducing integer valued variables l a , 

s - ^ E(^ - M + 2 - E n> + f E £(° 6 % - i E ^ + 2 - E E l % h % 

(A10) 

Now 6 a are not constrained and it is convenient to change variables to singlet and non-singlet 
components b° and b a , 

b a = fo t + b a f a (All) 

where t° and t a , a = 1..N — 1, are normalized generators of singlet and non-singlet symmetry 
groups respectively, i.e. £° = = and Y^aKfia = $ ab - We a ^ so introduce singlet 

and non-singlet combinations, 

/0 = ^E r ' r = E n « ( A12 ) 



Note that 1° and Z a are generally not integer valued. Hence, our action becomes, 

s - ^ E(^ - M + 2 ^°)l + f E( nfo % + 2 - E M - * E m 
+ |EE( Q& % + 2 -EE^ ( A13 ) 

a jfj, a jfj, 
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Shifting 27tV 'Nbf = dip + 2ir\/Nb } and integrating over ip, we obtain, 

+ |EE( n n>+ 2 -EE z fA> (A15) 

a jr> a 

with the constraint, 

(V ■ Z°)j = V^Vsj (A16) 



Moreover, performing a shift 6 a — > b a + in (|A14[) we learn that the flavoured components 
must satisfy, 

V ■ l a = (A17) 
The constraints (1A16I) and (1A17I) can be combined to give, 

(V-r)j = sj for all a (A18) 

Physically, l a corresponds to worldlines of vortices of e tda field. The singlet component of 
this field 1° corresponds to local vortices, when e l9a have the same winding for all a. Now, 
we can integrate the fields 6° and b a out in eq. (1A14I) . 



3 = 1^2*0^(2^ (A19) 

2 

+ ^E s ^'^+^E( n #^ ( A2 °) 

jj' " jj'n jn 

+ JEEWtf( 2T % ( A21 ) 



2 



with the kernels, 



u ii' 



D f 
u jj' 



_L \^ \ e iKj-j') CA22) 

V ^ e 2 N + gJ2 u 4sm 2 V ' 

I y 1 e *a-iO (A23) 

^V^E,4sin 2 ^f 1 J 



The first term in (IA19I) corresponds to a short range interaction between singlet components 
of vortices. This is in agreement with the expectation that local vortices should have local 
interactions. On the other hand, as we see from ( 1A21I) . non-singlet components of I, i.e. 
global vortices, have long-range interactions. 

Now let's see how our vortices couple to the magnetic source field Hj^. This coupling is 
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given by the second term in (1A19j) . Rewriting, 



= 2mJ2H ]lM (e 2 D^ f )(VNl ) I/tl = 2*1 £ H 3lt {^Df jt )(^ HjV (A24) 
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we observe that each vortex Z a carries a flux equal to 2ne 2 D s (k = 0) = 2?r/iV. This is again 
in accordance with our arguments in section IIHI 

Finally, for completeness we also discuss the terms in eq. (1A20I) . These terms are not 
very important as they don't couple to the dynamical variables l a . The first term in ( 1A20I) 



corresponds to bare, short range, interaction between monopoles, while the third term cor- 
responds to the monopole's contribution to the magnetic flux. 

Now, we would like to rewrite the action (1A19I) in a local form. We consider the action, 

s = i E E( D5 % + k E E(^ Q - Bat « - 2 i H+ 2 ™ Q )l {A25) 

»e(e^V; ( A26 ) 




The first line ( IA25j) is the simplest lattice generalization of the continuum action (13.91) . We 



have a set of A^ — 1 dynamical gauge field B a which couple to non-singlet combinations 
of currents associated with a set of A" matter fields V a ~ e %ipa . The integer variables 
ensure the periodicity of ip a . In the second line 0A260 we introduce the monopole source 
field s which couples to the combination Y2 a V 9 "; i n accordance with our continuum guess 
(13. 7p . Also, observe the coefficient 2tt/N, corresponding to the flux of a single vortex, in the 
coupling to the source field H. 

We have also introduced an additional term with coupling constant t in eq. (1A26j) . This 



term gives local interactions to the vortices of e tipa fields. We don't actually expect this term 
to drastically modify the critical properties of our theory (as we shall see, it will change some 
ultra-local interactions into just local ones), nevertheless, we have included it to make the 
similarity with the direct theory more pronounced. 

We would like to make the connection between actions (1A25|) and (lAlj) . As usual, we 



Poisson resum m a in (IA25I) with the help of integer valued variables I? and perform a shift, 



2-kw! = 2nm a - B a t a a + d(p a , 



s = i E E^b% + £f E !>"" - >>?, 4 E (e ' < A ^> 

~ 1 E ( E ¥?) s J + ^E l U W ° + 5 X - ^ a )j, (A28) 



By integrating over (p a we recover the constraint (lA18j) . Next, we go to the rotated variables 
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(1A12j) and perform analogous rotation on on m' a . Then, 



+ i E E< nB °)?, +^E E(™'°)i + < E E '!,( 2 ™'° + ( A3 °) 

a j/j, a jp, a 

Integrating over m' a we obtain, 

= (2tt£ y o _ 1^2 ^ ^ o ^ o (A } 

+ & E E(° B °)* + +i EE + 1 E E<ni (aso 

a jfj, a jfj, a 



We note that the above action, except for the last term in eq. (1A32|) . is exactly the same 



as that in (fATl) with the identification, m'° = b'°, B a = 2irb a , e 2 = g = 

t = g/N. As for the last term | Y^ a EE^i )? . ^ gives rise to an ultra-local interaction 
between global vortices. Since these vortices already interact with a long-range potential 
( 1A23I) . we don't expect this term to alter the critical properties of the theory. We also note 



that the somewhat unusual term with the coefficient t in (IA26I) is related to the coupling g 
in the sing let kernel D s flA22|) (but not in the flavoured kernel D s (TA"23~I) !). Setting t = 0, 
would make interactions between singlet vortices, ultra-local, instead of just short-range. 
Again, we don't expect such a change to alter the critical properties of the theory. 

Thus, we have argued that the lattice actions (1A25I) and ( jA If) are equivalent, up to local 
interactions in the flavoured sector. 
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